Excited state calculations using phaseless auxiliary-field quantum Monte Carlo: potential energy 

curves of low lying C2 singlet states 
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We show that the recently developed phaseless auxiliary-field quantum Monte Carlo (AFQMC) method can 
be used to study excited states, providing an alternative to standard quantum chemistry methods. The phaseless 
AFQMC approach, whose computational cost scales as M'^-M* with system size M, has been shown to be 
among the most accurate many-body methods in ground state calculations. For excited states, prevention of col- 
lapse into the ground state and control of the Fermion sign/phase problem are accomplished by the approximate 
phaseless constraint with a trial wave function. Using the challenging C2 molecule as a test case, we calculate 
the potential energy curves of the ground and two low-lying singlet excited states. The trial wave function is ob- 
tained by truncating complete active space wave functions, with no further optimization. The phaseless AFQMC 
results using a small basis set are in good agreement with exact full configuration interaction calculations, while 
those using large basis sets are in good agreement with experimental spectroscopic constants. 

PACS numbers: 71.15.-m,02.70.Ss,21.60.De,31.15.vn,31.50.Bc,31.50.Df,71.15.Qe 
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I. INTRODUCTION 

The ability to calculate electronic excited states of 
molecules and extended systems is necessary to predict key 
phenomena and properties of technologically important sys- 
tems. Compared to ground states, however, the accurate cal- 
culation of excitation energies is significantly more difficult. 
For molecules, a variety of many-body electronic structure 
quantum chemistry approaches have been developed, typi- 
cally using a one-particle basis to represent the many-body 
wave function. For small molecules with modest basis sets, 
the full configuration-interaction (FCI) method is exact, but 
FCI is not practical for realistic calculations, since the com- 
putational cost scales exponentially as the system size is in- 
creased. For larger systems, approximate coupled cluster 
(CC) methods' are the standard, but these methods also have 
rather steep computational scaling with system size {e.g., 
0{M'^) for CCSD(T), CC with single and double excitations 
and perturbational triplets, where M is the number of ba- 
sis functions]. For extended systems, less accurate approxi- 
mations based on density functional theory (DFT) and time- 
dependent DFT have been developed; GW and Bethe-Salpeter 
type methods have also been shown to be promising.^ Corre- 
lated quantum chemistry methods have also been embedded 
in DFT calculations to treat extended systems.-^- The most 
commonly applied quantum Monte Carlo (QMC) method in 
electronic structure has been diffusion Monte Carlo (DMC)r" 
which has also been used to compute excited statesi^ Com- 
pared to ground states, however, the accuracy of the results 
may depend on the symmetry^ and show greater sensitivity^ 
to the trial wave function used in the fixed-node approxima- 
tion to control the sign problem and maintain orthogonality. 

The recently developed phaseless auxiliary-field quantum 
Monte Carlo (AFQMC) metho d^i'"-"-'^ is an orbital-based al- 
ternative many-body approach. AFQMC can be expressed 
with respect to any chosen single-particle basis {e.g., gaus- 
sians, planewaves, Wannier, etc.), and it exhibits favorable 
©(M^-ikf) scaling. For ground states, the new AFQMC 
method has been applied to close to 100 systems, including 



first- and second-row moleculeS j"''^''-^ transition metal ox- 
ide moleculesjJS simple solidsrili post-rf elementsj'- van der 
Waals systemsr^ and in molecules in which bonds are being 
stretched or brokeniiii^ In these calculations we have oper- 
ated largely in an automated mode, inputting only the DFT 
or Hartree-Fock (HF) solutions as trial wave functions. The 
method demonstrated excellent accuracy, consistently able 
to correct errors in the mean-field trial wave function. In 
molecules, we have found that the accuracy of the phaseless 
AFQMC is comparable to CCSD(T) near equilibrium geom- 
etry and better when bonds are stretched. AFQMC thus pro- 
vides new opportunities for the efficient and accurate many- 
body calculations of ground and excited states in both molec- 
ular and extended systems. 

The seemingly simple C2 molecule presents a significant 
challenge for many-body methods.i^iaS The C2 molecule is 
difficult because of the strongly multireference nature of the 
ground state wave function [in which only 70% of the 
weight is the restricted Hartree-Fock (RHF) determinant] and 
the presence of nearby low-lying states. The shortcomings of 
standard quantum chemistry calculations for C2 were shown 
by recent benchmark FCI calculation&iS of the potential en- 
ergy curves (PECs) of its ^^S^ ground state and two low- 
lying singlet excited states. This benchmark shows that most 
correlated methods based on a single-determinant reference 
state wave function exhibit large nonparallelity errors 
(NPE — defined as the difference between the maximum and 
minimum deviations from FCI along the PEC). Spin-restricted 
CCSD(T) [referred to as RCCSD(T) hereafter] was found to 
exhibit a large NPE of 98 mi?h due to the poor behavior of 
RHF in the dissociation limit. Spin-unrestricted UCCSD(T), 
which is usually less accurate near equilibrium, has an NPE 
of 34m£'h. The excited state PECs are not accurately mod- 
eled by any of the commonly used single-reference methods, 
nor by CI including full quadruple substitutionSfi^ Similarly, 
a recent DMC study^ found that, even in its ground state at 
equilibrium geometry, the total energy of C2 showed a large 
fixed-node error ~ 40 mEh, if a single-determinant trial wave 
function is used. 
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As a new QMC method, the phaseless AFQMC provides an 
alternative route to the sign problem from fixed-node DMC. 
The random walks take place in a manifold of Slater determi- 
nants, in which fermion antisymmetry is automatically main- 
tained in each walker Applications have indicated that often 
this reduces the severity of the sign problem and, as a result, 
the phaseless approximation has weaker reliance on the trial 
wave function. It is interesting then to test the method for 
excited states, where QMC calculations depend more on the 
trial wave function and our knowledge of it is less. The chal- 
lenging Co molecule, where FCI results are available for the 
modest-sized basis set 6-3 IG*, provides an excellent test case. 

We first describe the AFQMC methodology for ground and 
excited states. We then make detailed comparions of our C2 
calculated results with the FCI calculations of Ref. |T^ Fi- 
nally, our calculated PECs and spectroscopic constants with 
large realistic basis sets are presented and compared with ex- 
perimental results. 



II. METHODOLOGY 

In this paper, we focus on the lowest-lying C2 singlet states: 
the X ^1;+ ground state and the B and B' ^Y.+ excited 
states. Since AFQMC uses a projection method to obtain the 
excited states, collapse to the ground state must be prevented. 
The B ^ A+ state belongs to a different irreducible representa- 
tion of the symmetry group of the Hamiltonian than does the 
X ground state, but the B' excited state belongs to 
the same irreducible representation as the ground state. Both 
of these cases are discussed below. We first briefly review the 
phaseless AFQMC method and then discuss the calculation of 
the excited states. 



A. Ground state 

Stochastic ground state quantum Monte Carlo (QMC) 
methodsri2i2L2^ which are exact in principle, use projection 
from any reference many -body wave function |$r), which 
has non-zero overlap with the ground state. In practice, how- 
ever, the Fermionic sign problem^i2i2i2i2^ must be controlled 
to eliminate exponential growth of the variance. For example, 
in DMC, a single- or multi-reference trial wave function is 
used to impose approximate nodal boundary conditions of the 
many-body wave function in electronic configuration space (a 
Jastrow factor is also included to reduce the stochastic vari- 
ance). By contrast, phaseless AFQMC samples the many- 
body wave function with random walkers { 1 0) } in the space of 
Slater determinants, which are expressed in terms of a chosen 
single-particle basis. Here we use standard quantum chem- 
istry gaussian basis sets.-S§^ Each |0) has the form of a HE or 
DFT wave function, with the orbitals varying stochastically 
in the projection. AFQMC controls the sign problem dif- 
ferently, using the complex overlap of the walker 10) with a 
trial/reference wave function 1$^) which is a determinant or a 
linear combination of determinants. 



The ground state energy is obtained from the mixed estima- 



tor 



En = 



(1) 



where H is the many -body Hamiltonian and I^Po) is the 
ground state wave function, which is given by imagi- 
nary time (/?) projection from Using the Hubbard- 
Stratonovich (HS) transformation^2i2i an importance sam- 
pling transformation- then expresses the mixed estimator as 
a stochastic average over the walkers and their Monte Carlo 
weights W0, 
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where the "local energy" El is defined as 



El[cI>] 



(2) 



(3) 



The energy computed from Eq. ^ is approximate. In ad- 
dition to a statistical error which can be accurately estimated 
and reduced with further sampling, there is a systematic error 
because of the phaseless constraint with |$r)- In other words, 
the realization of e~^^ in Eq. ([T]) is approximate. 



-f3H 



(4) 



because of the constraint on the random walk paths in the 
Slater determinant space (or equivalently, in the correspond- 
ing auxiliary-field space). This is the only approximation in 
the calculation. The computed ground state energy is not an 
upper boundi2i22 

In previous applications, phaseless AFQMC with a sin- 
gle unrestricted Hartree-Fock (UHF) determinant |<i>™^) was 
found to often give better overall and more uniform accu- 
racy than CCSD(T) in mapping PECs i"-'^''^ In some cases, 
however, such as the BH and N2 molecules, achieving quan- 
titative accuracy of a few mE^ for the entire PEC required 
multi-determinant We also find this to be true for the 

C2 molecule, as discussed below. Although a multi-reference 
containing Nd determinants increases the computational 
cost roughly by a factor of Nd, we have found that this is sig- 
nificantly offset by the gain in statistical and systematic accu- 
racy due to the use of a \^r) that more closely resembles the 
ground state 



B. Excited states 

Calculating the energy of the lowest excited state belong- 
ing to an irreducible representation that is different from that 
of the ground state [e.g., the B ^A^ state) is straightforward. 
In this case, we simply choose 1$^) in Eq- O to have the 
symmetry of the desired excited state. Since the Hamiltonian 
is invariant under the group of symmetry transformations, this 
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FIG. 1: (Color online) The overlaps of the AFQMC wave function with different states, and the computed total energy, as a function of 
imaginary time (5 at interatomic separation 1.25 A in C2. The left, center, and right panels show results for the X, B, and B' states, 
respectively, with AFQMC using the corresponding truncated CASSCF(8,16) wave function as |<l?r)- In each case, the upper panel shows 
the thi-ee overlap integrals \{$^ \'^mc)\^ , |("l?f |*mc)P, and IvI'mc)!^- The lower panels show the computed AFQMC energy, given 
by Eq. (|T](, together with the exact FCI energyi^ for comparison (indicated by the horizontal dashed lines). The |$r) are given by truncated 
CASSCF(8,I6) wave functions. The 6-31G* basis set is used. 



projects out, in Eq. Q, any component in the walker determi- 
nant 10) belonging to a different representation. Although it is 
usually not possible for a single-reference |$,.) to satisfy the 
symmetry requirement, multi -reference wave functions can, at 
least approximately. We use truncated complete active space 
self-consistent field (CASSCF)^ wave functions in this study, 
and our tests indicate that symmetry breaking due to the trun- 
cation is small (see below). 

For calculations of excited states belonging to the same ir- 
reducible representation as the ground state, e.g., the B' 
state, we rely on the fact that the corresponding refer- 
ence wave function is approximately orthogonal to the exact 
ground state, « 0. Obtaining accurate AFQMC 

results for excited states thus depends on using sufficiently ac- 
curate excited state trial wave functions. Our results indicate 
that a multi-reference with a modest number of determi- 
nants directly taken from a CASSCF calculation is adequate. 

Figure [T] illustrates our approach for the first three singlet 
states in C2 at a bond length near that of the ground state equi- 
librium. In AFQMC the wave function is given by 

where is the reference function used in the calculation for 
the phaseless constraint and for importance sampling, and the 
sum is over the random walker population, {w^, at each 
time slice. From Eq. (|5]), we can obtain estimates of the over- 
lap integrals | (<f>f,|^'Mc) P (where s ^ X, B, or B') to probe 



the composition of the AFQMC wave function. We normalize 
I'I'mc) by explicitly evaluating ^ (^'mcI^'mc) at each /?. Be- 
cause this involves "undoing" the importance sampling^' [di- 
vision by the factor ($r-|<?!') on the right-hand side of Eq. ©J, 
there are large statistical fluctuations, as can be seen in the up- 
per panels. An average population of 1000 walkers is used in 
these calculations. 

The left panel in Fig. [T] shows the ground state calculation, 
and its upper panel shows a large (~ 60%) | |^'mc}P 
overlap of the AFQMC walker population with the refer- 
ence wave function, as expected. By contrast, the overlap 
($f |5'mc)P is essentially zero, since the symmetry of the 
State is different from that of the \^^) ground state, 
and since symmetry breaking due to truncation of the full 
CASSCF wave funtion is evidently weak. Moreover, the 
somewhat larger I'I'mc)!^ overlap (~ 2 — 5%) is not 
surprising, since |$,^ ) has the same ^11+ symmetry as the 
ground state. The same trends are observed in the center and 
right panels of Fig. [T] In the center panel, where the total en- 
ergy of the B state is calculated, the overlaps of the walker 
population with the different symmetry X and B' states is 
extremely small. In each of the three panels, the walker popu- 
lation overlap with the corresponding trial reference state is 
large. Finally, it is interesting to note that the characteris- 
tic imaginary time (~ 20 — 40 E^^) to reach the asymtop- 
tic value of the overlaps is much larger than the equilibration 
times (~ 3E^^^) for the total energy. (Section [III Al presents 
detailed comparisons with the FCI results.) 
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FIG. 2: (Color online) Convergence of the AFQMC energy as a 
function of the number of determinants Nd in the multi-reference 
|$r), obtained by truncating a CASSCF wave function (see text). 
The upper panel shows the convergence for the X state near its equi- 
librium at i? = 1.25 A. The lower panel corresponds to the B state 
at i? = 1.8 A, a slightly larger bond length than the X-B crossover 
point at i? ~ 1.7 A). The arrows indicate the No given by our trun- 
cation criterion of 97% integrated weight. The 6-3 IG* basis set was 
used. Exact FCI energyl' is shown as horizontal dashed lines. 



C. Computational details 

Most of our AFQMC calculations use as |$,.) a truncated 
CASSCF(8,16) wave function, obtained from the GAMESS 
quantum chemistry programi^ The CASSCF wave function 
is truncated such that the weight (squared coefficient) of the 
retained determinants is ^ 97% of the total. Figure|2]plots the 
computed AFQMC energy as a function of the number Njj of 
retained determinants (ordered by decreasing weight) for the 
X and B states, calculated at i? = 1.25 and R = 1.8 A, re- 
spectively. Convergence to the 1 or 2 m£^h level is relatively 
quick, with the more correlated B state showing somewhat 
slower convergence. Our truncation criterion above gives 
Nd = 86 and Nd — 106 determinants for the X and B 
states, respectively. The corresponding AFQMC energies are 
seen to be well converged with respect to the truncation;^ 

Our AFQMC calculations use the local energy 
formalism^i^ using standard gaussian-type basis sets/^^ 
Reference wave functions were obtained using GAMESS^^ 
and the one- and two-body matrix elements were obtained 
using a modified NWCHEM code Ji A mean-field background 
subtraction is applied to the Hamiltonian prior to the HS 
transformation, which improves the computational efficiency 
and reduces systematic errorsjii^ In most of our AFQMC 
calculations, we use At = 0.01 E^^. We confirmed that the 
resulting Trotter error is less than 1 mEh by calculations at 
multiple At values at selected geometries in both the 6-3 IG* 
and larger basis sets. 

All runs use an average population of about 100 random 



walkers, with initial population generated from the RHF wave 
function or "broken symmetry" RHF»i^ (A short phaseless 
AFQMC projection is first invoked for P ^ 1 E^^, with 
the RHF wave function as \^r) and its copies to form the 
initial population; the resulting population, which is purely 
spin-singletji^ is then fed into the regular calculation with the 
CASSCF |$r) ) Typical runs have an equilibration phase of 
P ~ 10 E^^ and then a growth phase of (3 ^ 10 E^^^, in 
which the trial energy is adjusted via the growth estimator— 
and set for the rest of the simulation. A measurement of 
(3 ~ 150 E^^ is needed to achieve a statistical accuracy of 
1 mE^. To give an idea of the present computational cost, 
such a calculation with the cc-pVTZ basis at a single geometry 
(with Nd ~ 100 in the trial wave function) requires approxi- 
mately 2.5 days on one core of a 2.2GHz Opteron processor. 
The current implementation simply importsil one- and two- 
body gaussian matrix elements from other quantum chemistry 
programs. The Hamiltonian and overlap integrals are treated 
as dense matrices. No special properties of the underlying 
gaussian basis set are exploited. 

Our focus with the phaseless AFQMC has so far been on 
establishing the basic framework in a variety of systems, and 
testing its systematic accuracy and robustness. In this paper 
our main purpose is to present a proof of concept for excited 
state calculations. Although the present implementation (with 
gaussian basis sets and a density decomposition of the two- 
body interaction) has shown excellent accuracy, there remains 
considerable flexibility in the choice of the one-particle basis 
and the form of the HS transformation. It is possible that ex- 
ploiting the flexibility can lead to significant further improve- 
ment in accuracy and computational efficiency. 



in. RESULTS AND DISCUSSION 

As bonds are stretched in the dissociation of molecules, ac- 
curate treatment of strong electronic correlations becomes im- 
portant, especially in the intermediate regime when bonds first 
begin to break. The C2 molecule is a particularly challeng- 
ing example, and this trend is seen clearly below. The ground 
state at equilibrium geometry is already nontrivial»^ but as the 
bond is stretched in the ground state and across all geometries 
in the excited states, the systematic errors grow significantly in 
all calculations. We first compare our phaseless AFQMC re- 
sults with benchmark FCI calculations (More detailed dis- 
cussion about recent theoretical calculations on C2 molecule 
can be found in Refill) Calculated PECs and spectroscopic 
constants from realistic basis sets are then presented and com- 
pared with experimental results. 



A. Comparison with benchmark FCI results 

Figure[3]compares the phaseless AFQMC ground state PEC 
with the FCI results from Ref.[l9l using the 6-3 IG* gaussian 
basis set. In order to benchmark our AFQMC calculations, 
which do not employ the frozen core approximation, we es- 
timate a frozen-core correction to FCI using the difference 
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FIG. 3: (Color online) C2 ground state PEC; comparison of phase- 
less AFQMC using single-determinant UHF and two-determinant 
UHF2 reference wave functions (see text) with benchmark FCI re- 
sults from Ref.[l^. All calculations used the 6-3 IG* gaussian basis 
set. 

of UCCSD(T) energies with and without the frozen-core ap- 
proximation. The AFQMC/UHF PEC was calculated using a 
single determinant \^r) from a UHF calculation for the sin- 
glet ground state. The AFQMC/UHF2 PEC used the simplest 
multi-reference consisting of two determinants, 

IUHF2) cx |UHF) + a|UHFa,) , (6) 

where the parameter a is variationally optimized. In Eq. (|6]l, 
UHF) is the usual UHF wave function, which is a broken 
symmetry state with opposite spins on the two C atoms, but 
which preserves the axial symmetry of the molecule. The 
UHFai:} stated breaks the axial symmetry of the molecule 
and is analogous to the "antiferromagnetic solution" found in 
a local density approximation calculation in Ref. i39i. Near 
equilibrium, the energy of |UHF) is lower than that of 
UHFaa;), but, as the bond is stretched beyond i? ~ 1.5 A, the 
UHFjja;) energy becomes lower. The combined |UHF2) state 
in Eq. (|6]l has a variational energy that is ^ 10 mE\i lower than 
the I UHF) near equilibrium and ~ 20 mEh lower at i? ~ 1.8 

A. 

The AFQMC/UHF PEC has an NPE of 38m£;h while 
AFQMC/UHF2 has a smaller NPE of 20 mEh. The NPEs 
of standard quantum chemistry calculations are shown in Ta- 
ble m together with those of AFQMC. The single reference 
AFQMC/UHF NPE is thus seen to be significantly better than 
RCCSD(T) and comparable to UCCSD(T). Using only two 
determinants, AFQMC/UHF2 has a smaller NPE than the 
CISDTQ result. (AFQMC/UHF2 has a lower energy at large 
bond lengths due to the change in the leading terms of the FCI 
wave function from the single RHF-like configuration near 
equilibrium to a two-determinant configuration in the dissoci- 
ation limit.) Nevertheless, an NPE of ~ 20mi?h is unaccept- 
ably large for a high-level method such as QMC. We also note 



TABLE I: Nonparallelity error (NPE) of standard quantum chem- 
istry methods for the C2 ground state PEC (taken from Ref.fT^. com- 
pared with that of the phaseless AFQMC method using UHF, two- 
determinant UHF2, and truncated CASSCF(8,16) trial wave func- 
tions. The range of bond lengths is i? = 0.9 - 3.0 A. All calculations 
used the 6-3 IG* gaussian basis set. 



Method 


INrn (mil/h) 


RHF 


212 


UHF 


78 


MP2 


130 


RCCSD(T) 


98 


UCCSD(T) 


34 


CISD 


116 


CISDT 


51 


CISDTQ 


26 


AFQMC/UHF 


38 


AFQMC/UHF2 


20 


AFQMC/CASSCF 
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that removing spin-contamination in the walker population, as 
discussed in Ref. [T^ does not yield significant improvements 
to either AFQMC/UHF or AFQMC/UHF2. This indicates that 
these I <!>,,) are themselves poor. 

We now show that the phaseless AFQMC results become 
much more accurate when multi-reference |$r) are used. As 
described in Sec. Ill CI these are truncated CASSCF(8,16) 
wave functions, with no further optimization. Figure |4] com- 
pares the AFQMC/CASSCF and FCI calculated PECs. We 
see that all three PECs are mapped out accurately, including 
the B' PEC which is of the same symmetry as the ground 
state X. The overall accuracy of the AFQMC PECs for both 
the ground state and the excited states is better than Bmi^h 
for all but one point, the smallest bond length in B. The X 
and B crossover at i? ~ 1.7 A and the B and B' crossover at 
i? ~ 1.1 A are both accurately described. With the truncated 
CASSCF trial wave function, the energies also appear to be 
variational in all cases, while the AFQMC/UHF energies are 
not. 

Tablellllpresents the spectroscopic constants corresponding 
to the PECs in Fig.H] Since CCSD(T) is poor in the dissoci- 
ation limit, the dissociation energies in Table ilD are calcu- 
lated using the energy of the free C atom for each method. For 
AFQMC, this somewhat improves the comparison of with 
FCI, since the AFQMC energy of the C atom is more accurate 
than that from the value in Fig. |4] in the dissociation limit, 
where the truncated CASSCF(8,16) |$r) is not very good. 
The spectroscopic constants obtained from energy expectation 
values of the truncated CASSCF |<i>r) have substantial fitting 
uncertainties, which originate from noise in the determinant 
truncation. Both the AFQMC/CASSCF and the full CASSCF 
results are in very good agreement with the exact FCI results, 
while those based on energy expectation values of the trun- 
cated CASSCF \^r) are significantly worse. This shows that 
AFQMC improves significantly over the truncated CASSCF 
trial wave functions. 
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*(C-C) (Angstrom) 

FIG. 4: (Color online) AFQMC/CASSCF PECs for the three 
lowest lying singlet states in C2, using multi-reference truncated 
CASSCF(8,16) I*,.), compared with FCI results.J^ FCI results are 
given by the solid curves and AFQMC results are given by sym- 
bols with enoi bars. RCCSD(T) results for the ground state are also 
shown as the orange line. The inset shows deviations from FCI in 
mEh- All calculations used the 6-3 IG* gaussian basis set. 



TABLE II: AFQMC/CASSCF calculated C2 spectroscopic con- 
stants, corresponding to Fig. |4] using multi-reference truncated 
CASSCF(8,16) For comparison, results from FCI,— CCSD(T), 
Ihe full CASSCF, and from the energy expectation values of the trun- 
cated CASSCF |<I>,.) are also shown. All calculations used the 6- 
31G* basis set. Results for the equilibrium bond length (in A), 
vibrational frequency uoe (in cm"'^), and ground state dissocation en- 
ergy De (in eV) are shown. For excited states, the excitation energy 
Te is defined as the energy difference between the minima of the 
excited and ground states. Combined statistical and fitting errors in 
AFQMC are shown in parantheses, while pure fitting uncertainties 
are shown in square brackets. 

CASSCF(8,16) 





full 


truncated CCSD(T) 


QMC 


FCI 




ground state 










re 


1.2563 


1. 27 1 [4] 


1.2577 


1.2567(3) 


1.2581 




1893 


1705 [76] 


1869 


1888(12) 


1863 


De 


6.530 


5.12[1] 


5.953 


6.085(5) 


6.030 


B'At 


excited state 












1.3997 


1.39[3] 




1.4004(9) 


1.4000 




14I4[1] 


I360[332] 




1394(13) 


1391 


Te 


1.712 


1.75[8] 




1.874(6) 


1.761 




excited state 










Te 


1.3868 


I.399[7] 




1.4009(12) 


1.3931 




1463 


1453 [127] 




1325(32) 


1394 


Te 


1.920 


1.91[3] 




2.189(6) 


2.058 



B. Realistic C2 results using large basis sets 

In this subsection, we present phaseless AFQMC/CASSCF 
PECs with large basis sets and compare with experimen- 
tal spectroscopic values. Figure |5] shows AFQMC/CASSCF 
PECs for the three lowest lying singlet states, using truncated 
CASSCF(8,16) wave function as in the cc-pVTZ basis 




*(C-C) (Angstrom) 



FIG. 5 : (Color online) C2 AFQMC/CASSCF PECs using | $^ ) given 
by a truncated CASSCF(8,16) wave function. For comparison, the 
RCCSD(T) ground state PEC is also shown. All calculations used 
the cc-pVTZ basis.— 



setis£ For comparison, the figure also shows the RCCSD(T) 
ground state PEC, which was calculated using the NWCHEM^l 
computer program. Table Ulllpresents the corresponding spec- 
troscopic constants and experimental values. We have also 
computed the ground state spectroscopic constants using the 
larger cc-pVQZ basis set^^ shown in Table |IV| As seen from 
these two tables, the residual finite basis set error is small and 
the computed spectroscopic constants are nearly converged. 
(Ground-state results obtained using the cc-pV5Z basis set by 
the contracted multi-reference CI (CMRCI) method^ are also 
included in Table llVI ) 

The calculated AFQMC results in Tables |III] and |IV] are 
in very good agreement with experimental values. The 
CASSCF(8,16) results, i.e., using the full CASSCF WF are 
also quite good. Each CASSCF(8,16) WF consists of 414864 
determinants in the D2h space group. The truncated CASSCF 
used in AFQMC, however, uses only the first 52 - 275 
determinants, which amounts to 97% of the total CASSCF 
determinant weight. At smaller distances, the correlation ef- 
fects are small and thus fewer determinants are needed. At 
large bond lengths, there are more determinants included in 
the |$r), indicative of stronger correlation effects. As seen 
in Table HIH AFQMC significantly improves the results when 
compared to those obtained from the variational estimate us- 
ing the truncated CASSCF WF. The level of agreement here 
between the AFQMC results with larger basis sets and exper- 
iment is consistent with that in the smaller basis set between 
AFQMC and FCI in the previous subsection. 



IV. SUMMARY 

We have shown that molecular excited-state calculations 
are possible with the phaseless auxiliary-field QMC method. 
Using CASSCF trial WFs, the method delivers accurate PECs 
in the challenging C2 . The computed spectroscopic constants 
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TABLE III: AFQMC/CASSCF calculated C2 spectroscopic con- 
stants compared to experiment. Results from CCSD(T), the full 
CASSCF, and from the energy expectation values of the truncated 
CASSCF wave functions are also shown. Conventions are as in Ta- 
ble |ll] All calculations used the cc-pVTZ basis set. 



for the lowest three singlet states, two of which have the same 
symmetry, are in very good agreement with experiment. 



CASSCF(8,16) 

full truncated CCSD(T) 



X ^ ground state 

re 1.2479 

UJe 1862 

De 6.53 



QMC Expt. 



1.262[3] 
1766[67] 
4.64[1] 



1.2508 
1842 
6.03 



B '^A^ excited state 

r-e 1.397[1] 1.416[7] 

LOe 1395[23] 1351[101] 

Te 1.511[5] 1.38[3] 



1.2462(9) 1.2425 
1884(17) 1855 
6.32(1) 6.33 



1.391(1) 1.3855 
1376(23) 1407 
1.723(7) 1.498 



B 



' ^Eg excited state 
r-e 1.381 



1.398[4] 
LOe 1489[4] 1392[63] 

Te 1.779[2] 1.68[2] 



1.393(1) 1.3774 
1441(12) 1424 
2.082(8) 1.910 



TABLE IV: AFQMC/CASSCF calculated C2 ground state spectro- 
scopic constants compared to experiment. Conventions are as in 
Table Ull Calculations used the cc-pVQZ basis set (except CMRCI 
which used the cc-pV5Z basis). 

CASSCF(8,16) 
full truncated CCSD(T) CMRCI" QMC Expt. 



X Eg ground state 

Te 1.2452 1.262[3] 1.2459 

uje 1868 1759[29] 1852 

De 6.57 4.69[1] 6.19 



1.2467 1.244(1) 1.2425 
1853 1850(21) 1855 
6.29 6.41(1) 6.33 



"Values from analytical fitting in Ref.l4Q 
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